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ABSTRACT 

G54. 1+0.3 is a young pulsar wind nebula (PWN), closely resembling the Crab, for which no thermal 
shell emission has been detected in X-rays. Recent Spitzer observations revealed an infrared (IR) 
shell containing a dozen point sources arranged in a ring-like structure, previously proposed to be 
young stellar objects. An extended knot of emission located in the NW part of the shell appears 
to be aligned with the pulsar's X-ray jet, suggesting a possible interaction with the shell material. 
Surprisingly, the IRS spectrum of the knot resembles the spectrum of freshly formed dust in Cas A, 
and is dominated by an unidentified dust emission feature at 21 /im. The spectra of the shell also 
contain various emission lines and show that some are significantly broadened, suggesting that they 
originate in rapidly expanding supernova (SN) ejecta. We present the first evidence that the PWN 
is driving shocks into expanding SN ejecta and we propose an alternative explanation for the origin 
of the IR emission in which the shell is composed entirely of SN ejecta. In this scenario, the freshly 
formed SN dust is being heated by early-type stars belonging to a cluster in which the SN exploded. 
Simple dust models show that this interpretation can give rise to the observed shell emission and the 
IR point sources. 

Subject headings: 



1. INTRODUCTION 

G54. 1+0.3 is a pulsar wind nebula (PWN) with prop- 
erties very similar to the Crab. Radio observations show 
a nebula with a flat sp ectru m, approximately 2' in ex- 
tent dReich et al.lll984t lGreedll985t iVelusamv fc Becker! 
1988), and no evidence of the supernova remnant (SNR) 
shell. While the extent of the X-ray nebula in the Crab 
is 3 times smaller than in the radio, non-thermal X-ray 
emis sion from G54. 1+0.1 extends out to the radio bound- 
ary ((Seward 1989; Lu et alj |2001). Chandra observations 
of G54. 1+0.3 revealed a point source surrounded by a 
ring and torus, and a nebula elongated in the E/W direc- 
tion (jLu et al.ll2002D . All components have a power-law 
spectrum with an index steepening with increasing dis- 
tance from the pulsar. There is no evidence for a thermal 
component associated with supernova (SN) ejecta or the 
swept up interstellar medium (ISM). 

The 137 ms pulsar J1930+1852 that pow ers the neb- 
ula was discovered by iCamilo et alj (|2002f ) , who calcu- 
late a characteristic age of 2900 yr. Considering a pos- 
sible range o f braki ng indices and initial spin periods, 
ICamilo et all (|2002D estimate the age of G54. 1+0.3 to 
be between 1500 and 60 00 yr. From the equations de- 
scribing PWN evolution, iChevalierl (|2005l ) calculates an 
age of 1500 yr and and an initial spin period Pq of 100 
ms. Based on HI line emission and absorption measure- 
ments, the distance to G54.1+0.3 is in th e 5-9 kpc range 
(jWeisberg et al.l[200l[Leahv et aUl2008h . while the pul- 
sar's dispersion measure implies a d istance less than or 
equal to 8 kpc (|Camilo et al.ll2002ft . iLeahv et al.l (j2008ft 
recently suggest a morphological association between the 
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nebula and a CO molecular cloud at a distance of 6.2+0.1 
kpc. However, the absence of X-ray thermal emission and 
the lack of evidence for an interactin of the SNR with 
the cloud argue against this association. In this paper, 
we assume a distance of 6 kpc for G54. 1+0.3, at which 
10"«0.3pc. 

While there has been no previous evidence for an inter- 
action of the PWN with SN ejecta, circumstellar material 
(CSM), or the ISM in G54.1+0.3, a recent Spitzer dis- 
covery of an infrared (IR) shell surrounding the PWN 
may be the first evide nce of such an interaction (Slane 
120081 : IKoo et al.l 12008ft . The first IR detection from 
G54. 1+0.3 was made with the Infrared Astronomical 
Satellite (IRAS) at 25 fim, 60 /mi, and 100 /tm, and 
the emission was consistent wi th a cold dust component 
with a temperature of ~ 30 K (jArend t 1989; Sakc n et alJ 
11992ft . Spitzer imaging revealed an IR shell surrounding 
the PWN in which a dozen point sources are arranged 
in a ring-like structure, suggested to be young stellar ob- 
jects (YSOs) whose formatio n was triggered in the late 
stages of the progenitor's life (|Koo et al 1 12008) . The mor- 
phological association between G54. 1+0.3 and the shell 
is evident in the IR images and suggests that the PWN 
may be interacting with the shell material. The nature 
of the shell still remains unclear. While the morphology 
suggests that the shell is SN ejecta and freshly formed 
dust swept up by the PWN, presence of YSOs would 
require that it be a preexisting circumstellar shell or a 
molecular cloud. In this paper, we present deep Chan- 
dra imaging of G54. 1+0.3 and Spitzer IR imaging and 
spectroscopy of the associated shell that provide new ev- 
idence that the shell emission originates from SN ejecta 
and dust. 

The paper is organized as follows; Section [2] describes 
the observations and reduction of the Chandra and 
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TABLE 1 

CHANDRA AND SPITZER OBSERVATIONS 

47.9s 



Detector/ Prog. Observation Exposure 
Channel ID Date Time 



Chandra 

ACIS-S 9886 2008 Jul 8 66.17 ks 

9108 2008 Jul 10 35.11 ks 

9109 2008 Jul 12 164.33 ks 





9887 


2008 Jul 15 


25.16 ks 


Spitzer 


IRAC 


3647 


2005 Oct 21 


2 frames X 30 s 


MIPS 24 


3647 


2005 May 15 


20 cycles X 30 s 


IRS SL 


40736 


2007 Nov 05 


5 cycles X 14 s 


IRS LL 


40736 


2007 Nov 05 


5 cycles X 6 s 


IRS SH 


40736 


2007 Nov 05 


10 cycles x 120 s 


IRS LH 


40736 


2007 Nov 05 


5 cycles X 6 s 



Note. — The MIPS 70 )im image was obtained from 
the MIPSGAL survey ICarev et alj|200SI) . 



Spitzer data, Section [3] describes the basic analysis and 
fitting of the X-ray spectra and IR line emission, Section 
IHpresents a detailed analysis of the IR spectroscopy, Sec- 
tion [5] discusses the dust emission, Section |6] describes a 
new interpretation for the origin of the IR shell and the 
embedded point sources, and finally, Section [7] summa- 
rizes the conclusions of our study. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Chandra 

Observations of G54. 1+0.3 were taken with the Ad- 
vanced CCD Imaging Spectrometer, ACIS-S, on board 
the Chandra X-ray observatory on 2008, July 8, 10, 12, 
and 15, under the observation ID numbers 9886, 9108, 
9109, and 9887. The observations were carried out using 
the FAINT mode and the corresponding exposure times 
were 66.17, 35.11, 164.33, and 25.16 ks, for a total expo- 
sure time of 290.77 ks. The standard data reduction and 
cleaning were performed using Ciao Version 3.4. 

The cleaned Chandra datasets were merged into a sin- 
gle event file using the merge_all task in Ciao version 
3.4. The final merged image is shown in Figure [TJ Spec- 
tra were extracted from the merged file from eight dif- 
ferent regions covering the PWN and one background 
region. The ex traction regions are similar to those in 
ILu et all (|2002D for the purpose of comparison, and are 
shown in Figure [2] Because the observations were car- 
ried out close in time and have the same chip orientations 
and positions on the sky, the corresponding effective area 
and spectral response files were generated from a single 
observation with ID 9109. Source spectra were back- 
ground subtracted and grouped to include a minimum of 
75 counts in each bin for region 1, 50 counts in each bin 
for regions 2, 3, 4, 7, and 8, and 20 counts in each bin 
for regions 5 and 6. The observations are summarized in 
Table [U 

2.2. Spitzer 

Infrared spectroscopy of the IR shell was carried out 
with all modules of the Infrared Spectrograph (IRS) 
(jHouck et all 12004ft aboard Spitzer on 2007, Nov. 05, 
covering the 5-37 /jm wavelength range (program ID 
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Fig. 1.— Chandra ACIS image of G54. 1+0.3 in the 0.3-10.0 keV 
band, shown on a linear scale (0 to 120 counts). 

40736). The location of the shell with respect to the 
PWN is shown in Figure [31 and the positions of the IRS 
slits are shown in Figure |U The low resolution spectra 
were taken at position 19h30m28.0s, +18°52'15.60", cen- 
tered on the diffuse shell emission, but also covering the 
bright knot in the northwest with the long-low (LL) mod- 
ule. The low resolution observations had 5 cycles for each 
module, with an exposure time of 14 s and 6 s each for the 
short-low (SL) and LL modules, respectively. The high 
resolution spectra were taken at two different positions; 
the diffuse shell emission at 19h30m28.0s, +18°52'15.60", 
and the bright knot at 19h30m26.40s, +18°52'07.00". 
The locations of the slits for the background spectrum 
are also shown in Fig |U The short-high (SH) spectra 
were obtained using 10 cycles of 120 s, and the long-high 
(LH) using 5 cycles of 6 s each. 

Data were processed with pipeline version S17.0.4 and 
cleaned using IRS Rogue Pixel Mask Editing and Im- 
age Cleaning software (IRSCLEAN1.9). The background 
was subtracted from the high resolution data, individual 
Basic Calibrated Data (BCDs) of each nod were median- 
combined, and spectra were extracted with the Spitzer 
IRS Custom Extractor (SPICE) using full slit extrac- 
tions with the extended source calibration for both low 
and high resolution modules. The subsequent analysis, 
including averaging of the individual nods, background 
subtraction for the low resolution modules, trimming of 
the order edges, and fitting of the emission lines, was 
carried out using the Spectrosc opic Modeling, An alysis 
and Reduction Tool (SMART) (jHigdon et al.112001 . For 
the low resolution modules, the spectra extracted from 
the off-source nods were used for the background. We 
have also included the imaging data from the Infrared Ar- 
ray Camera (IRAC) and Multiband Imaging Photometer 
(MIPS) aboard Spitzer, shown in Figures [3] and [5j The 
IRAC observations at 3.6, 4.5, 5.8, and 8.0 /im were car- 
ried out on 2005, Oct. 21, under the program ID 3647, 
using two 30 s frames. The data were processed with the 
pipeline version S14.0.0 and reduced with the MOsaicer 
and Point source Extractor (MOPEX) version 18.1.5. 
The MIPS 24 /im observation was carried out on 2005, 
May 15, under the program ID 3647, using 20 cycles of 
30 s, and wa s processed with the pipeline version SI 6. 1.0 
(|Slandl2008ft . The MIPS 70 fim image of G54. 1+0.3 was 
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Fig. 2. — Extraction regions for the X-ray spectra, overlaid on the 
smoothed Chandra ACIS image of G54.1+0.3 in the 0.3-10.0 keV 
band (shown on a 0—230 logarithmic scale). The fitted parameters 
for each of the numbered regions are summarized in Tabled 



obtained from the MIPSGAL survey (|Carev et alJl2009h . 
Spitzer observations are summarized in Table [TJ 



3. ANALYSIS 



3.1. 



X-ray 

The deep Chandra ACIS image of G54. 1+0.3 in the 
0.3-10 keV band is shown in Figure [TJ The central pul- 
sar is surrounded by a 5'.' 7 by 3'.'7 X-ray ring, with the 
long axis oriented roughly north-so uth, suggesting an in- 
clination angle of about 40 degrees (|Lu et alj r2002). The 
pulsar and the ring are embedded in a diffuse nebula, 
2'.0 by 1 f 3 in size, which appears to have a bipolar elon- 
gation, running roughly east-west, perpendicular to the 
apparent plane of the ring. 

Spectra were extracted from the Chandra data using 
the apertures s hown in Figure El The regions are sim- 
ilar to those in iLu et a l. (2002) and include the pulsar 
in region 1, the surrounding ring in region 2, the west- 
ern arc in region 3, the western jet in region 4, two 
bright knots of emission in the east in regions 5 and 
6, and fainter extended emission between these struc- 
tures and in the outer regions of the PWN in regions 
7 and 8, respectively. The background spectrum was 
extracted from an elliptical region with an area of 7710 
square arcseconds, centered at the position 19h30m29.6s, 
+ 18°50'35.25". The spectrum for each of the regions 
was fitted with an absorbed power law model using the 
Ciao 3.4 Sherpa software. The fitting of region 1 also 
included a pileup model, since the point source count 
rate of 0.077 s _1 is expected to caus e significant pileup 
that affects the spectral parameters (Davis 2001). The 
absorbing column density was fit simultaneously for re- 
gions 2-8, and the column density for the region encom- 
passing the point source was then fixed to this value. 
The fitting results are summarized in Table El and the 
spectra are shown in Figures [6] and [71 We find an absorp- 
tion column density of (1.9 5 ± 0.04) x 10 22 cm -2 , some- 
what greater than that of ILu et al.l (|2002f ), who found 
Nh = (1.6 + 0.1) x 10 2 2 cm -2 . We find a similar value to 
that of lLu et alj (|2002l ) if we include region 1 in the joint 
fit and do not account for pileup. When the spectrum 
of each region is fit separately, the Nh (10 22 cm -2 ) in re- 
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Fig. 3.— Three color image of G54.1+0.3, where the MIPS 24 
fim emission is shown in blue, MIPS 70 fim emission in green, and 
radio emission in red. The Chandra X-ray contours are shown in 
white. 

gions 2-8 ranges from 1.88-2.19, while the best fit for 
region 1 gives an JVjy(10 22 cm~ 2 ) of 1.54. This suggests 
that a more accurate estimate for Nh is achieved if the 
piled-up region 1 is excluded from the fit. 

The spectral fits show that the outer regions of the dif- 
fuse emission in G54. 1+0.3 have a softer spectrum than 
the interior regions. We find no evidence of a thermal 
component associated with shock-heated ejecta or ISM 
in the X-ray spectra. The new estimate of the absorb- 
ing column density results in a slightly higher photon 
index for the pulsar, T = 1.44 + 0.04. The ring and E/W 
elongations have similar photon indices, which are harder 
(smaller) than those in the outer diffuse nebula, suggest- 
ing lower synchrotron losses or more particle injection 
in these regions. There is no evidence for emission lines 
in any of the spectra. The parameters derived from the 
initial analysis of the Chandra data are used in our anal- 
ysis of the IR observations. A more detailed study of the 
X-ray data will be presented in a separate publication. 

3.2. IR Morphology 

The Spitzer IRAC and MIPS images of G54. 1+0.3 re- 
veal an infrared shell south and west of the PWN, with 
the X-ray nebula filling the cavity of the shell. The shell 
is detected at 5.8, 8.0, 24, and 70 fim, and is approxi- 
mately 1 '5 in radius at 70 fim. A shell structure is evi- 
dent in the 5.8 /im and 8.0 /urn IRAC images, shown in 
Figures [5ji and [SJd, along with significant emission that 
extends far from the PWN to the north that may be as- 
sociated with Polycyclic Aromatic Hydrocarbon (PAH) 
emission unrelated to G54. 1+0.3. The structure of the 
shell in the MIPS images appears to have a a brighter 
western lobe, and a fainter southeastern lobe (Figure [5J; 
and[3Ji). The X-ray and radio emission from the PWN 
fills the cavity of the shell. The western jet structure, in 
particular, appears to extend directly into a void region 
in the IR nebula, terminating at the brightest region in 
the 24 fim MIPS image, i.e. the IR knot. This is best 
shown in Figure [3j where the MIPS 24 /im emission is 
shown in blue, MIPS 70 fim emission in green, radio 



4 



Temim et al. 



TABLE 2 

CHANDRA ACIS SPECTRAL FITTING RESULTS 



Region Description Area(arcsec 2 ) Photon Index ATjy(10 22 cm 2 ) Fx (observed) Fx (unabsorbed) Reduced \ 2 



1 


Pulsar 


22.03 


1.44 ± 0.04 


(1.95) 


2.10E-12 


3.26E-12 


1.03 


2 


Ring 


90.77 


1.86 ± 0.05 


1.95 ± 0.04 


5.88E-13 


1.21E-12 


0.66 


3 


West arc 


230.9 


1.89 ± 0.05 




6.64E-13 


1.40E-12 


0.92 


■1 


West jet 


185.4 


1.79 ± 0.06 




4.06E-13 


7.88E-13 


1.08 


5 


East knot 1 


16.94 


1.90 ± 0.12 




8.94E-14 


1.90E-13 


0.63 


6 


East knot 2 


19.12 


1.97 ± 0.13 




7.62E-14 


1.73E-13 


0.70 


7 


Inner nebula 


1586 


2.05 ± 0.04 




1.84E-12 


4.45E-12 


0.94 


8 


Outer nebula 


5559 


2.20 ± 0.04 




1.18E-12 


3.34E-12 


0.83 



Note. — The extraction regions are shown in Figure [2] Listed uncertainties are 1.6 cr (90 % confidence) statistical uncer- 
tainties from the fit. The Nh value listed for Region 2 is the result of the simultaneous fit for regions 2—8. The fluxes were 
calculated in the 0.3-10 keV band and are in the units of erg cm" 2 s _1 . 



TABLE 3 
IRAC AND MIPS FLUXES 



Wavelength 


Flux 


Extinction Corrected 


W 


(Jy) 


Flux (Jy) 


5.8 


~ 0.2 


~ 0.3 


8.0 


~ 0.6 


~ 1.0 


24 


23.7 ± 2.4 


40 ± 4 


24 - IR knot 


3.8 ± 0.4 


6.5 ± 0.7 


70 


76 ± 15 





Note. — IRAC fluxes were estimated based on 
the average flux density at selected positions scaled 
to the size of the shell. The uncertainties include 
IRAC and MIPS calibration uncertainties, but do 
not account for the uncertainties in the extinc- 
tion correction. The extinction co rrection was ap- 
plied using the extinction curve of Chiar & Tielensl 
d2006T) . 

emission in red, and Chandra X-ray contours in white. 
The most interesting features of the 24 /im image are 
11 point sources embedded in th e diffuse shell th at are 
arranged in a ring-like structure (|Koo et al.l l2008). The 
positions of the point sources are marked in Figure 3J 
The two northern point sources are outside of the ring 
structure, but appear to be located at the tips of faint 
ridges of emission that trace back to the shell structure. 
The IRAC and MIPS fluxes from the IR shell and the 
knot are listed in Table [3] Since the shell region in the 
IRAC images is crowded with stars, the fluxes were es- 
timated by scaling the average flux density from several 
positions on the shell by the total area of 4 square ar- 
cminutes. In measuring the flux from the IR knot, we 
used a background region from the diffuse shell emission 
east of the knot. 

3.3. IR Spectroscopy 

The high and low resolution IR spectra of G54. 1+0.3 
are shown in Figures [8] and [9j and the positions of the 
IRS slits are shown in Figure HI The low resolution spec- 
trum shows a rising continuum longward of 15 /im, and a 
broad emission feature peaking at approximately 21 /im. 
This feature is more pronounced in the high resolution 
spectrum of the bright knot, while it is less pronounced 
in the spectrum from the diffuse shell (Figure E|) . The 
high resolution spectra also reveal a broad emission fea- 
ture around 12 /im, and a narrow emission feature around 
11.3 /im that may be associated with PAH emission or 
other carbonaceous grains. The high resolution slit from 
which the background spectrum was extracted is located 




Fig. 4. — Positions of the IRS slits overlaid on the MIPS 24 /im 
image. The LL slits are the wide slits oriented in the East/West 
direction. The SL slits are oriented North/South. The LH and 
SH slits are the smaller slits in yellow, positioned at three different 
locations; the diffuse shell emission at the same location as the 
low resolutions slits ( position 1), the extended IR knot in the west 
(position 2), and the background region located off of the shell to 
the west. The point sources embedded in the diffuse shell that are 
evident in the MIPS 24 /an image are marked by blue circles. 



in a region of low background emission, and it underesti- 
mates the contribution from the emission in the vicinity 
of the IR shell that is evident in the IRAC 8 /im image 
(FigurelSJa) . However, the background regions for the low 
resolution spectrum are more representative of the local 
background emission. Since the low resolution spectra 
were extracted from the entire SL and LL slits, they are 
not characteristic of the point-to-point spectrum along 
the slits that varies spatially. The sharp peak around 7.5 
/im is an artifact caused by a mismatch between orders, 
and was trimmed from the residual source spectrum. 

There are a total of ten emission lines present in the 
spectrum of the IR shell; the molecular hydrogen H2 S(2) 
(12.3 ^m), and H 2 S(l) (17.0 /im) lines, and ionic lines 
of [Ar II] (6.99 /im), [S IV] (10.51 /im), [Ne II] (12.81 
Aim), [CI II (14.37 /im), [S III] (18.71 /im), [Fell] (25.99 
/xm), [S III] (33.48 /im), and [Si II] (34.82 /im). The po- 
sition of the line at 14.37 /im that we identify as [CI II] 
is near the position of [Ne V] 14.32 /im but since we 
do not detect Ne III in our spectrum, we are confident 
that the line is produced by chlorine. The line-fitting pa- 
rameters for the LH and SH modules are listed in Table 
21 where positions 1 and 2 correspond to spectra from 



"G54.1+0.3 with Chandra, Spitzer 



5 



TABLE 4 
IRS LINE FITS 



Line ID 


Line Center 


Line Flux 


De-reddened Flux 


FWHM 


FWHM 


True FWHM 


Shift 




w 


(10 -6 erg /cm? /s/sr) 


(10~ 6 erg/cm 2 /s/sr) 


w 


(km/s) 


(km/s) 


(km/s) 



Position 1 

H 2 (12.2786) 
[Ne II] (12.8135) 
[CI II] (14.3678) 
H 2 (17.0348) 
S III] (18.7130) 
Fe II] (25.9883) 
S III] (33.4810) 
Si II] (34.8152) 

Position 2 

S IV] (10.5105) 
Ne II] (12.8135) 
CI II] (14.3678) 
H 2 (17.0348) 
[S III] (18.7130) 
[S III] (33.4810) 
[Si II] (34.8152) 



12 


2775 


± 





0013 


1.24 ± 0.31 


2.15 ± 0.54 





018 


± 





001 


435 ± 24 




-29 ± 35 


12 


8138 


± 





0008 


11.15 ± 0.76 


18.72 ± 1.28 





024 


± 





001 


568 ± 24 


279 ± 104 


5 ± 20 


14 


3721 


± 





0005 


4.59 ± 0.16 


7.79 ± 0.27 





049 


± 





001 


1029 ± 21 


902 ± 32 


89 ± 9 


17 


0385 


± 





0004 


2.74 ± 0.10 


5.10 ± 0.19 





030 


± 





001 


525 ± 18 


177 ± 165 


65 ± 6 


18 


7195 


± 





0010 


13.51 ± 0.83 


26.61 ± 1.64 





051 


± 





001 


823 ± 16 


658 ± 44 


104 ± 16 


26 


0048 


± 





0084 


4.62 ± 1.26 


7.33 ± 1.99 





095 


± 





008 


1089 ± 92 


966 ± 33 


191 ± 97 


33 


4890 


± 





0026 


25.28 ± 2.67 


35.94 ± 3.79 





075 


± 





003 


675 ± 27 


450 ± 70 


72 ± 23 


34 


8174 


± 





0040 


72.69 ± 6.67 


102.76 ± 9.43 





136 


± 





004 


1172 ± 34 


1059 ± 30 


19 ± 35 


10 


5164 


± 





0019 


1.07 ± 0.19 


2.47 ± 0.43 





021 


± 





001 


596 ± 28 


334 ± 87 


168 ± 55 


12 


8123 


± 





0003 


9.74 ± 0.30 


16.36 ± 0.50 





029 


± 





001 


675 ± 23 


450 ± 70 


-28 ± 6 


14 


3689 


± 





0011 


3.12 ± 0.28 


5.28 ± 0.48 





037 


± 





001 


780 ± 21 


604 ± 48 


23 ± 21 


17 


0392 


± 





0016 


1.18 ± 0.22 


2.20 ± 0.40 





028 


± 





002 


490 ± 35 




77 ± 29 


18 


7212 


± 





0006 


33.92 ± 1.49 


66.81 ± 2.94 





041 


± 





001 


651 ± 16 


424 ± 69 


132 ± 9 


33 


4939 


± 





0011 


48.08 ± 2.68 


68.36 ± 3.80 





059 


± 





001 


528 ± 9 


160 ± 198 


116 ± 10 


3d 


8229 


± 





0046 


57.30 ± 7.56 


81.00 ± 10.68 





108 


± 





005 


934 ± 43 


786 ± 40 


66 ± 40 



Note. — Listed uncertainties are l-cr statistical uncertainties from the fit and do not include the IRS wavelength and flux calibration 
uncertainties. 



TABLE 5 

IRS LINE FITS IN LL & SL OVERLAP REGION 



Line ID 



Line Center (/im) Line Flux De-reddened 



Region 1 

Ar II] (6.9853) 
Ne II] (12.8135) 



CI II 
S III 
S III 



(14.3678) 
(18.7130) 
(33.4810) 



Si II] (34.8152) 



6.993 ± 0.001 
12.813 ± 0.013 
14.401 ± 0.013 
18.733 ± 0.004 
33.539 ± 0.029 
34.889 ± 0.015 



84.9 ± 4.9 
6.4 ± 1.9 
3.2 ± 1.3 
11.9 ± 1.1 
33.7 ± 8.4 
79.3 ± 10.6 



122.4 ± 7.1 

10.8 ± 3.3 
5.4 ± 2.2 
23.5 ± 2.2 

47.9 ± 11.9 
112.1 ± 14.9 



Note. — Fluxes are in units of 10 — 6 erg/cm 2 /s/sr. Listed un- 
certainties are l-tr statistical uncertainties from the fit and do not 
include the IRS wavelength and flux calibration uncertainties. 



the diffuse sheil and the bright IR knot, respectivefy (see 
Figured]). Line-fitting parameters from the overlapping 
region of SL and LL modules are listed in Table El The 
extinction co rrection was app li ed us ing the IR extinc- 
tion curve of iChiar fe Tielensl (|2006l ) , and the relation 
N H /A K = 1.821 xlO 22 cm" 2 (|Drainelll989D . As shown in 
Figure El background emission significantly contributes 
to the observed lines of Ne, S and Si, implying that the 
uncertainties on their measured line intensities are prob- 
ably larger than the statistical uncertainties quoted in 
Table S) A number of emission lines in the high resolu- 
tion spectrum show evidence of broadening with respect 
to the spectral resolution of IRS of A/ A A ~ 600, includ- 
ing CI, S, Fe, and Si. The broadened line profiles of the 
Si and S lines are shown in Figure ITD1 

For the purpose of examining the spatial variation 
in the emission line intensities across the LL slits, 
we created 2-dimensional spectral line maps using the 
CUbe Builder for IRS S pectral Maps (CUBISM) soft- 
ware (|Smith et all [2001 . The maps were made by se- 
lecting a desired wavelength region for the source flux 
and subtracting a wavelength-weighted background, av- 
eraged across a specified wavelength region. The maps 



were created by selecting a source flux from 20.9-21.9 
/im, with a continuum from 20.5-20.7 /mi and 22.3-22.6 
/jm to show the spatial variation of the 21 /im feature; a 
source flux between 18.4-18.9 /im with a continuum be- 
tween 18.0-18.4 /im and 18.9-19.0 /tm for the [S III] line 
at 18.7 /tm; a source flux between 33.1-33.9 /im, with a 
continuum between 31.9-32.9 /im and 33.9-34.2 /tm for 
the [S III] line at 33.5 /im; and a source flux from 34.6- 
35.2 /tm, with a continuum from 34.1-34.4 /tm for the 
[Si II] line at 34.8 /tm. The narrow wavelength ranges 
for the backgrounds were chosen to avoid other emis- 
sion features and large gradients in the continuum slope. 
We smoothed the spectral line maps by a 2-pixel Gaus- 
sian before plotting the intensity profiles as a function 
of position along the slit (Figure fTTj). Spectra were also 
extracted from 7 spatial positions along the LL slit in 
3x2 pixel sub-slits. Due to the low extraction area, only 
the brightest lines of S and Si were measurable in these 
spectra. Their intensities are listed in Table [6j where the 
position along the LL slit is denoted by the pixel number 
on the x-axis in Figure 1111 

4. LINE EMISSION 

Spectra of the diffuse shell were obtained with all four 
modules of the IRS (Figure 2]). The high resolution SH 
and LH spectra were taken at two different positions on 
the shell; at the position of diffuse emission near the in- 
terface between the shell and the PWN (position 1), and 
at the position of the bright IR knot in the northwest 
(position 2). The LL slit was oriented in the East/ West 
direction and it spans the entire length of the shell, in- 
cluding the IR knot. The SL slits were roughly perpen- 
dicular to the LL slits, and span the length of the shell 
in the North/South direction. 

4.1. Line Profiles and Broadening 

The high-resolution spectra at both positions in the 
IR shell show evidence of spectral line broadening. The 
fitted line parameters for the two positions are summa- 
rized in Table |4j The expected resolution of LH and 
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Fig. 5. — Spitzer imaging of the shell surrounding the PWN in G54.1+0.3. The shell is detected in the IRAC 5.8 /im (panel a) and 8.0 
jim (panel b) bands. The MIPS 24 fim and 70 fim images are shown in panels c and d, respectively. The IR knot is indicated by the white 
arrow. The overlaid contours represent the Chandra ACIS contour of the PWN. 



SH is A/AA i~ 600, equal to a full width at half maxi- 
mum (FWHM) of 500 kms" 1 . The [CI II], [Fe II], and 
[Si II] lines at position 1 show the largest broadening 
of over 1000 kms -1 . The [S III] lines at this position 
are broadened to approximately 700-800 kms -1 . The 
[Ne II] line at position 2 shows evidence for slight broad- 
ening, while the widths of the H 2 lines are approximately 
equal to the resolution of the IRS. It is likely that in- 
adequate subtraction of the local background from the 
LH and SH spectra contributes significantly to spectral 
lines with narrow widths. The true FWHM in Table 2] 
was calculated by subtracting the instrumental FWHM 
of the IRS from th e observed FWHM in quadrature (see 
iDasvra et al]|200g ). 

The lines in the spectrum at the position of the IR 
knot (position 2) are broadened by a lesser amount than 



at position 1. The [CI II] and [Si II] lines again show 
the greatest broadening of ~ 800-900 km s -1 , while the 
[SIII] lines are only slightly broadened to an average of 
~ 600 km s -1 . We also note that the best-fit centroid of 
the [S III] lines at both positions appears to be redshiftcd 
by an average of ~ 120 kms -1 . Although this shift is 
not very large considering that the uncertainty of the IRS 
wavelength calibration for the high-resolution modules is 
on the order of 100 km s -1 , it is consistently larger than 
the shift of all the other observed lines. The exception 
is the [Fe II] line with a redshift of 191 ± 97 kms -1 , 
but this line is faint and the shift has a high uncertainty. 
The profiles of the lines with the highest broadening show 
evidence for structure. The high resolution line profiles 
of [Si II] 34.8 /im and [S III] 33.5 fim at both positions are 
shown in Figure ITUl for comparison. At position 1, the 
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Fig. 6. — Chandra X-ray spectra of the PWN in G54. 1+0.3 extracted from regions 1-4, shown in Figure|2] The spectra in all regions are 
well fit by an absorbed power-law model with the parameters listed in Table [2] 



[Si II] line appears to have a brighter component that 
peaks at approximately the same wavelength as [S III], 
in addition to a fainter blue shifted wing. 

4.2. Spatial Variation in Line Intensities 

The variation in the surface brightness across the LL 
slit of the [S III] 18.7 /tm, [S III] 33.5 fim, and the [Si II] 
34.8 /im lines, and the SL spectral line map of the [Ar II] 
6.99 /im line are shown in Figure [TT1 The peak surface 
brightness for each line has been normalized to unity, so 
their relative fluxes are not represented by the plot. The 
intensities as a function of position for all lines in Figure 
ITTI show structure that correlates with the MIPS 24 /tm 
image of the shell. The intensity of the [S III] 33.5 /tm 
line is well correlated with the 24 /tm emission and it 
appears to stay roughly constant across the western part 



of the IR shell, with a possible slight rise at the location 
of the bright IR knot. The spatial intensity profile of the 
[S III] 18.7 /im line is similar to the profile of the 21 /tm 
feature and shows a sharp peak at the position of the IR 
knot, implying an increase in density at this position (see 
Section fQ)) . 

Curiously, the intensity of the [Si II] line has a large 
peak at the apparent interface between the PWN and 
the IR shell. This peak is a factor of two brighter than 
the [Si II] emission across the rest of the shell and does 
not correlate with spatial features shown by the IR con- 
tinuum emission or with the [S III] lines. The relative 
intensities of the [S III] and [Si II] across the LL slit are 
also summarized in Table [6j where the line fluxes were 
measured at 7 positions from spectra extracted from 3x2 
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Fig. 7. — Chandra X-ray spectra of the PWN in G54.1+0.3 extracted from regions 5-8, shown in Figure|2] The spectra in all regions are 
well fit by an absorbed power-law model with the parameters listed in Table [2] 



pixel sub-slits. The highest line broadening is seen at 
roughly the same position where the [Si II] line peaks, 
position 1 of the high-resolution spectrum. As shown in 
Figure [TTJ the intensity of [Ar II] as a function of posi- 
tion along the SL slit also peaks near the same position. 
Figure shows the spatial profiles of the 24 jum shell 
emission and the [Ar II] and [Nc II] lines across the SL 
slit. The intensities of 24 /im emission and [Ar II] appear 
to be anti-correlated spatially. The [Ar II] and [Ne II] 
lines also appear to be anti-correlated across the SL slit. 
The Ne line is brightest along the northern part of the 
IR shell, positions 0-13 in Figure [T2J Figure [T3] shows 
sample spectra from three different positions across the 
SL slit, extracted from 2 by 6 pixel sub-slits. The three 
positions show spectra with weak [Ne II] emission, strong 



[Ne II] emission, and emission with prominent broad fea- 
tures at 9 and 12 pm. The spectrum with the strong 
broad features was extracted from the position of an em- 
bedded point source, suggesting that these features may 
be associated with the source, or enhanced in its vicinity. 



4.3. Evidence for Multiple Emission Regions 
The IRS spectrum shows that the [CI II], [Fe II] 



and 



[Si II] lines arc all broadened by a similar amount at both 
positions of the high-resolution slits. Their FWHM, de- 
convolved by the instrumental resolution (true FWHM 
in Table 2} , is on the order of 1000 km/s. Another simi- 
larity is that these lines all have low ionization potentials 
of less than 13.6 eV. The [S III] lines show less broaden- 
ing, have centroids that are slightly redshifted relative to 
the other observed emission lines, and have a higher ion- 
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TABLE 6 
IRS LINE FITS ACROSS LL SLIT 





Line ID 


Line Center 


Line Flux 


De-reddened 


SIII Line Ratio 


Electron Density 






(cm) 


(10 -6 erg/cm? / s/ sr) 


(10 -6 erg/cm 2 / s/ sr) 


(18.7/^m/33.5/im) 


(cm -3 ) 


Region 1 (Pixel 


14) 












[Si II] (34.8152) 


Q c; ni _|_ 


n no 


35.4 ± 4.5 


49.9 ± 6.4 






Region 2 (Pixel 


17) 












[S III] (18.7130) 


18.76 ± 


0.01 


4./ ± U.o 


y.o it i.o 






S III (33.4810) 


QQ R1 _|_ 

oo.Di it 


U.UZ 


9.5 ± 1.6 


13.5 ± 2.3 


0.69 ± 0.16 


< 800 




Si II] (34.8152) 


qa n^ J_ 


n no 
U.UZ 


47.6 ± 7.2 


67.7 ± 10.1 






Region 3 (Pixel 


20) 












[S III] (18.7130) 


18.73 ± 


0.01 


Q 1 i. Q 1 

y.i ± z.l 


ion 4. ^ q 

lo.U it 4.o 






[S III] (33.4810) 


QQ C'7 _|_ 


n no 


15.4 ± 3.8 


21.8 ± 5.4 


0.83 ± 0.28 


< 1300 


[Si II] (34.8152) 


q a nn _l 
o4.yu ± 


n m 


106.3 ± 8.4 


149.9 ± 11.9 






Region 4 (Pixel 


23) 












[S III] (18.7130) 


18.73 ± 


0.01 


i o i -i- i n 

1Z.1 it l.U 


oq q _i_ o n 
zo.o it z.U 






[S III] (33.4810) 


QQ ro _|_ 


n no 


28.9 ± 6.7 


41.0 ib 9.5 


0.58 ± 0.14 


< 600 


[Si II] (34.8152) 


1A QO 4- 


n m 


78.2 ± 10.7 


110.3 ± 15.1 






Region 5 (Pixel 


26) 














S III] (18.7130) 


18.73 ± 


0.01 


ofi c; _i_ o n 
zu.O it Z.\j 


c;o i _i_ /i n 
Oz.l it 4.U 








S III] (33.4810) 


33.44 ± 


0.02 


29.6 ± 8.1 


42.0 ± 11.5 


1.24 ± 0.35 


500 - 2400 




Si II] (34.8152) 


34.89 ± 


0.02 


64.6 ± 8.6 


91.7 ± 12.3 






Region 6 (Pixel 


29) 














S III] (18.7130) 


18.72 ± 


0.01 


23.9 ± 3.1 


47.0 ± 6.0 








S III] (33.4810) 


33.55 ± 


0.01 


54.9 ± 7.7 


77.9 ± 10.9 


0.60 ± 0.11 


< 600 




Si II] (34.8152) 


34.84 ± 


0.01 


50.5 ± 7.4 


71.6 ± 10.6 






Region 7 (Pixel 


32) 












[S III] (18.7130) 


18.71 ± 


0.01 


11.1 ± 3.7 


21.8 ± 7.3 






[S III] (33.4810) 


33.53 ± 


0.02 


21.1 ± 4.2 


30.0 ± 6.0 


0.73 ± 0.28 


< 1100 


[Si II] (34.8152) 


34.90 ± 


0.03 


48.5 ± 7.8 


68.8 ± 11.1 







Note. — Listed uncertainties are l-cr statistical uncertainties from the fit and do not include the IRS wavelength and 
flux calibration uncertainties. The ranges of estimated densities assume a temperature range of 100-10000 K, and take 
into account the statistical uncertainties on the [S III] line ratios. 



ization potential of 23.34 eV. The intensities of the low 
ionization lines do not appear to correlate spatially with 
the intensities of S and Ne lines. While the Si emission 
strongly peaks at the position where the PWN appears 
to encounter the IR shell, the spatial plot of the S and 
Ne intensities shows no enhancement at this position. 
We are not able to determine if the S and Ne lines are 
spatially correlated, since the SL and LL slits sample 
different regions of the shell. The spatial distribution of 
the Ar line is anti-correlated with Ne, and it appears to 
peak at approximately the same region where Si is en- 
hanced. It may be that the spectral line emission from 
the IR shell arises from distinct regions; a region where 
the highly broadened lines of ions with low ionization po- 
tentials are produced, including Si, CI, Fe and perhaps 
Ar, a region from which most of the sulfur emission orig- 
inates, and a region giving rise to the [Ne II] line that 
may or may not be spatially correl ated with sulfur. Sim - 
ilar components were observed bv lNeufeld et al.l (2007), 
who analyzed IRS spectral line maps from several SNRs 
interacting with molecular clouds, and found that the 
line emission generally originates from spatially distinct 
groups. The groups included lines from ions with low ion- 
ization potentials such as [Fe II] , [Si II] , and [P II] , lines of 
ions with ionization potentials higher than 13.6 eV, not 
including [S III], and the lines of [S III]. While the lines 
in various groups sometimes showed similar spatial dis- 



tributions, the lines within a gr oup were always strongly 
correlated (jNeufeld et al. 1 12007} ). The emitting regions in 
G54. 1+0.3 may represent distinct spatial components of 
different temperatures, or arise due to different critical 
densities of the ionized species. A possible explanation 
for the origin of different emission regions will be dis- 
cussed in Section l6.2.1l 

4.4. Electron Density 

The electron density of the emitting material can be 
probed using the temperature-insensitive density diag- 
nostic ratio of the ground state fine structure lines of 
[S III] at 18.7 fun and 33.5 /im. Both lines are present in 
the LL module, and the ratio of their surface brightness 
maps is shown by the dashed curve in Figure 1111 Since 
the spectral maps were used to calculate the [S III] line 
ratios across the slit, the absolute values of the point-to- 
point ratios in Figure [TT1 are highly uncertain. Since the 
line profiles were not fitted in this case, the shape of the 
underlying continuum and the variation of the [S III] line 
strength in the local background all contribute to the un- 
certainty. However, it is clear that the ratio peaks at the 
location of the bright IR knot, implying that the elec- 
tron density here sharply increases. We estimated the 
electron densities using the extinction corrected intensi- 
ties of the fitted [S III] lines in the LH and SH spectra 
(Tableffl and at multiple positions across the LL slit (Ta- 
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Fig. 8. — High resolution Spitzer IRS SH and LH spectra of the IR shell in G54. 1+0.3. The fainter spectrum is from position 1, the 
interface between the PWN and the shell, and the brighter spectrum is from position 2 at the IR knot (see Figure [4). 



ble ■ The calculation wa s performed using the IRAF 
STSDAS Nebular package (IShaw fc Dufouri fl995h . The 
high resolution [S III] 18.7 ^m/33.5 fim line ratios are 0.7 
± 0.1 and 1.0 ± 0.1 for positions 1 and 2, respectively, 
giving a density range of < 750 cm -3 at position 1, and 
500-1300 cm -3 at position 2, for an assumed tempera- 
ture range of 1000-10000 K. The quoted uncertainties are 
statistical uncertainties from the fit only. The ratios have 
additional uncertainties due to the fact that the LH and 
SH slits have different sizes and orientations, which af- 
fects the average surface brightness in the slit. Since the 
local background also contains [S III] lines, spatial vari- 
ations in the background add additional unquantifiablc 
uncertainties to the line ratios and the density estimates. 
Table [5] shows the [S III] ratio at six different positions 
across the LL slit. The ratios for all regions are consis- 
tent with being at the low density limit, except for the 
region corresponding to the IR knot at which the density 
is obviously higher, in the 500-2400 cm~ 3 range. 



4.5. Shock Diagnostics 

The PWN in G54. 1+0.3 appears to fill the cavity of the 
IR shell, with the pulsar jet terminating at the IR knot 
(Figure [3]). Due to a morphological correlation between 
the PWN and the IR shell, it is possible and even likely 
that the pulsar wind is sweeping through the shell mate- 
rial. A similar interaction is observed in the Crab Neb- 
ula, in which the PWN drives shocks into the SN ejecta 
and produces optic al emission from a shel l at the outer 
edge of the PWN (|Sankrit k. Hesterlll997l ). In order to 
estimate the velocity of the shock propagating into the 
shell material, we compared the measured spectral line 
intensities with predictions from models of shock waves 
in SN ejecta. Our model code does not include chlorine. 
However, CI and Ar should be found in the same region 
in the ejecta, and the excitation rates of the lines are sim- 
ilar, so the ratio is mainly determined by the abundance 
ratio. In addition, some emission may also originate in 
unshocked ejecta, and these ejecta may emit more effi- 



G54.1+0.3 with Chandra, Spitzer' 



11 



10 
8 

>; 6 



4 

2 




i 1 1 — i 1 — i — i 1 1 1 1 i — i 1 — i 1 1 r 



[Sill 





I I I 



15 



20 25 30 

Wavelength (/U-m) 



35 




8 10 
Wavelength (/i.m) 



12 



14 



Fig. 9. — Low resolution Spitzer IRS spectra extracted from the entire LL slit (top panel) and SL slit (bottom panel). The solid grey 
spectra represent spectra before background subtraction, dashed grey spectra are from the background extracted from the off-source nod 
positions, and the black spectra are the background-subtracted source spectra. The peak at around 7.5 /im is due to an order mismatch at 
the edges, and was trimmed in the final spectrum. 



ciently in [CI II] than in [Ar II] (see Section 16.2. 1|) , as 
Cl + is more likely than Ar + to be the dominant ion- 
ization species in freely-expanding unshocked ejecta. In 
the overlap region (Table [5]), the [Ar II] line is 23 times 
as bright as [CI II] , implying an abundance ratio similar 
to lWooslev fe W eaver (1995) nucleosynthesis models for 
core-collapse SNe. We use the [Ar II] /[CI II] ratio in the 
low-resolution to estimate the [Ar II] intensities at Posi- 
tions 1 and 2 of the high-resolution spectrum and scale 
all the intensities to [Si II] = 100 for ease of comparison 
in Table [3 

Models of emission from sho cks in SN e jecta were 
run using the c ode d escribed in iBlair et aLl (|2000l ) and 
iWilliams et al.l (|2008l ). The results are summarized in 
Table [71 These models differ from models of shocks with 
normal astrophysical abundances in that the enormous 
cooling rates generally give T e much lower than Ti in 
most of the cooling zone. The models assume the abun- 



dance set H : He : C : N : O : Ne : Mg : Si : S : 
Ar : Ca : Fe = 12.0 : 10.9 : 10.0 : 10.0 : 14.0 : 14.2 : 
14.3 : 14.3 : 14.0 : 14.2 : 13.0 : 13.3, given on a logarith- 
mic scale, but only the abundances of Ne, Si, S, Ar and 
Fe can really be constrained by the observations. The 
pre-shock density is taken to be 30 cm~ 3 , as suggested 
by the measured density of the IR knot. The density in 
the rest of the shell is somewhat lower, but this should 
only affect the density sensitive [S III] ratio. The gas is 
followed until it cools to 300 K, and the pre-shock ioniza- 
tion state is taken to be in equilibrium with the radiation 
field produced by the shock. 

The models should not be considered unique, in 
that the pre-shock ionization state, the abundances of 
coolants that are not constrained by the observations (es- 
pecially O), and the low temperature cutoff of the models 
will all affect the predicted spectra. Moreover, it is un- 
likely that a single shock speed really accounts for all the 
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TABLE 7 

SHOCK MODEL LINE INTENSITIES 



Shock Speed Ne II Ne III Si II 

(km/s) 



S III S III 

(18 fim) (33 fim) 



S IV Ar II Ar III Fe II 



Model Intensities 

20 8 
25 24 
30 42 



100 
9 100 
101 100 



Observed Intensities 

Position 1 18 < 2 

Position 2 20 < 2 



100 
100 



5 
47 
194 



26 
83 



5 
43 2 
150 32 



35 
81 



< 2 
3 



11 
80 
118 



86 
74 



1 5 
18 8 
115 11 



< 2 

< 2 



7 

< 2 



1.0 



0.8 



Note. — Shock m odels were run using the code described in Blair ct al. (200Cj) and 
IWilliams et~ al. (2008) and abundaces H : He : C : N : O : Ne : Mg : Si : S : Ar : Ca : Fe 
= 12.0 : 10.9 : 10.0 : 10.0 : 14.0 : 14.2 : 14.3 : 14.3 : 14.0 : 14.2 : 13.0 : 13.3. For comparison 
with models, we normalized the observed intensities to [Si II] 35 fim = 100. The [Ar II] line 
is not measured at Positions 1 or 2, so we assume that the [Ar II] /[CI II] ratio is constant. 
The preshock density of the models is 30 cm -3 . 
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Fig. 10. — High resolution spectral profiles of the [Si II] and [S III] 
33.5 /im lines at two positions in the IR shell. The 1-cr statistical 
uncertainties on each data point are on the order of 20%. 

emission. Therefore, we can draw only the rough con- 
clusions that the intensity ratios (or limits on ratios) of 
[Ne III]/[Ne II], [S IV] /[S III] and [Ar III]/[Ar II] suggest 
a shock speed of about 25 km s -1 or a little less. The 
abundances of Ne, Si, S and Ar are comparable, while 
that of Fe is about 10 times smaller. We do not really 
have a good handle on the O abundance, because the 
shock speed is too low to produce much [O IV] . Overall, 
the emission lines indicate a shock moving through the 
ejecta of a fairly massive SN, though a more O-rich com- 
position cannot be excluded. The high pre-shock den- 
sity in the IR knot, and possibly in the rest of the shell, 
suggests that shocks are driven into denser-than-average 
ejecta clumps (see Section RT21 . similar to what has been 
inferred for SNR B0540-69!3 lWilliams et alll200l ). 

5. DUST EMISSION 

IR imaging and spectroscopy of G54. 1+0.3 can provide 
important information about the properties and distri- 
bution of dust in the IR shell. The study of ISM dust 
provides constraints on models describing dust composi- 
tion and properties. These models are best described by 
a grain composition that predominantly consists of as- 
tronomical silicates and carbonaceous material, such as 



PAHs, graph ite, and amorphous carbon (for a review see 
lDraindl2009f h While a large fractio n of interstellar dust 
is grown m the ISM fD raind [2009h . the rest is injected 
by stellar winds, PNe, novae, and SNe. Dust is formed 
in the cooling gas, where the density is high enough for 
grain growth, and the dust composition depends on the 
available gas abundances and environmental conditions 
in these objects. Models show that SN ejecta can have 
a range of compositions, with some of the most abun- 
dant grains being those of MgSiOs, S1O2, Mg2SiQ4, Si, 
and C (for a review see lKozasa et al.ll2009t l. The shape 
of the IR spectrum and the presence of various emission 
features are clues that can help determine the mass and 
composition of dust in G54. 1+0.3. 

5.1. The Unidentified 21 /im Feature 

The IR spectrum of the shell (Figure[5]) shows a sharply 
rising continuum that flattens out beyond 21 /im. The 
most obvious difference in the high resolution spectra 
at the two positions in the shell is the strength of the 
broad emission feature at 21 /jm. The feature is much 
more pronounced in the IR knot than in the rest of the 
shell. This is best seen in Figure [TT1 where the relative 
intensity of the emission in the 20.9-21.9 /jm range is 
shown as a function of the position along the LL slit. 
The precise shape of the spatial intensity plot is affected 
by the shape of the underlying continuum, but it is ob- 
vious that the peak intensity of the 21 /jm feature is an 
order of magnitude brighter at the IR knot. The feature 
peaks at approximately 20.8 /im and has a FWHM of 
~ 3.7 /im, although this estimate is sensitive to the de- 
termination of the underlying continuum emission. The 
shape is asymmetrical; it rises steeply at shorter wave- 
lengths, and falls off more gradually longward of 20.8 /im. 
Additional features observed in the IR spectra that may 
or may not be associated with the IR shell are emission 
features at ~ 8.6 /jm and ~ 12.4 /im, with a FWHM of 
approximately 1.4 /im and 2.1 /i,m, respectively. These 
features are evident in Figure [H] and Figure 1131 

The 21 /im feature in the IR shell of G54.1+0.3 is re- 
markably similar to th e emission featur e observed in the 
IR spectrum of Cas A (|Rho et al.ll2008h . When we over- 
lay our high-resolution spectrum from position 1 with the 
IRS spectrum of Cas A shown in Figure 3 of iRho et alj 
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Fig. 11. — Intensity of various emission features in the G54. 1+0.3 spectrum as a function of the position across the low-resolution IRS 
slits. The position of the LL slit, represented by the blue box, is overlaid on the MIPS 24 /an image of G54. 1+0.3. The pixel scale is 
5.1"/pfxel. The curves in the plot show the spatial intensity variation of the 21 fim feature (green), and the [Si II] (blue), [S III] (red & 
orange) lines. The peak intensity of each feature has been normalized to unity and the amplitudes do not represent actual flux values. The 
black dotted line shows the actual ratio of the fluxes of the [S III] 18.7 /an and 33.5 fim lines. The overlaid slit running N/S represents 
a normalized surface brightness map of the [Ar II] line at 6.99 fim. The statistical 1-cr uncertainties of the intensities at each pixel are 
approximately 0.05 for the 21 /an feature, 0.09 for [Si II], 0.09 for [Ar II], and 0.12 for [S III] lines. 




(2008), wc find that the features peak at the same wave- 
length and have the same width, implying that they 
likely arise from the same dust composition. In addition, 
the spectrum of Cas A also shows the broad emission 
features at 9 fim and 12 /im that are also evident in our 
spectrum, suggesting that these emission features may be 
correlated with the 21 fim feature a nd produced by the 
same grain species. iRho et all ([2008) attributed the dust 
that produces the 21 fim feature to freshly formed dust 
in the SN ejecta, and as we will discuss in Section [6l we 
believe that the IR emission in G54. 1+0.3 also originates 
from freshly-formed SN dust. 

The IR spectrum of Cas A was fitted with multiple 
grain compositions, with the main contribution from 
Si0 2 , Mg protosilicates, and FeO (jRho et all I2009D . 
While the overall continuum was well fitted by this com- 
position, the fit had significant residuals from the 21 fim 
feature and from the entire broad emission feature at 12 
fim. The Si02 grains modeled as a continuous distribu- 
tion of ellipsoids (CDE) provided a somewhat better fit 
to the 21 fim feature in Cas A, in addition to being able to 
produce t he 9 fim feature when modeled at a higher tem- 
perature (|Rho et al.ll2009T l . In our attempts to fit the 21 
fim feature in the spectrum of G54. 1+0.3, we found that 
the peak produced by the Si02 CDE composition was 
too narrow. However, if the Si02 grains are porous, the 
feature may be broadened and shifted to better match 



the observed profile. For this reason, these grains still 
may be plausible candidates for the production of the 21 
fim feature in Cas A and G54. 1+0.3. 

In order to gain insight into other possible grain species 
that could give rise to the 21 fim emission, we turned to 
the extensive work that has been done in identifying the 
"21 fim feature" in carbon-rich protoplanet ary nebulae 
PPNe). The emission was discovered by iKwok et al.l 
1989) from IRAS observations of PPNe, and has since 
been observed in other dusty environments, including 
planetary nebulae with Wolf Rayet p rogenitors and tw o 
extreme carbon stars (for a review see lPosch et al.ll2004j) . 
The emission feature observed in these objects actually 
peaks at a wavelength of 20.1 fim, somewhat lower than 
the feature in G54. 1+0.3. In a recent Spitzer study of 
carbon-rich PPNe, high resolution spectroscopy revealed 
that the 21 fim emission is accompanied by a narrow 
emis sion feature at 11.3 f im and a broader feature at 12.3 
fim (jHrivnak et al.|[2~009f ). similar to the feature in Cas A 
and what we observe in the IR shell. This provides some 
evidence that the features observed at various bands are 
related and that perhaps the same species give rise to 
emission in PPNe and in the IR shell of G54. 1+0.3. 

While various grain species have been proposed as 
the carriers of the 21 fim resonance feature, including 
TiC, SiS 2 , Si0 2 , SiC, and FeO, all but the latter two 
have been discarded based on abundance constraints and 
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the presence of other unobserved resonance b ands (e.g., 
iZhang et all 120091: ISpeck fc Hofmeister] |2004[ ) . In a re- 
cent paper, IZhang et all (|2009f ) ruled out all the men- 
tioned carrier candidates, except FeO, due to an insuffi- 
cient abundance of Ti, S, or Si in sixteen PPNe. They 
found that the FeO nano dust can reproduce the 21 /im 
feature in these sources without producing unobserved 
secondary features nor exceeding the Fe budget. While 
FeO is a good fit to the 21 /im feature in PPNe, the peak 
wavelength of 19.9 /jm falls short of the 20.8 /jm peak ob- 
served in G54. 1+0.3 and Cas A. When modeled with FeO 
grains with a CDE distribution, the peak wavelength of 
the feature shifts to 20.7 //m, but the FWHM increases to 
6.5 /im, far too broad to fit the feature we observe. FeO 
also fails to produce the other observed emission features 
at 9 /im and 12 /im that are likely related to the 21 /im 

fe ature. 

ISpeck k. Hofmeisten ()2004f ) suggested that the 21 /im 
feature may be produced by SiC grains that contain 
impurities, and through laboratory experiments showed 
that doped SiC grains do indeed produce a resonance at 
21 /im that matches the shape of the feature in PPNe. 
They also found that nano-SiC grains produce a peak 
that is shifted redward by ~ 1 /im, matching the peak 
of the featu re observed i n the IR shell of G54. 1+0.3 
and Cas A. ISpeck et al.l (|2005| ) found that amorphous 
and nanocrystalline samples all produce three peak com- 
plexes near 9, 11, and 21 /im. The SiC grains seem ca- 
pable of producing all dust emission features that we 
observe in the shell of G54. 1+0.3, including the sharp 
emission feature at 11.3 /im, the broader emission fea- 
tures around 9 and 12 /im, and the 21 /im feature with 
the same peak wavelength as we observe. The relative 
strength of these features depends on the grain size and 
shape, and the temperature of the underlying contin- 
uum. The strength of the 21 /im feature in particular, 
likely depends on the amount of carbon impurities in the 
SiC sample. 

Since strong emission from the 21 /im feature is 
mostly obse r ved in the PPN phase of stellar evolution, 
ISpeck et al.l {2005) suggest a scenario in which the ap- 
pearance of this feature is due to the particular physical 
conditions in these objects, rather than a production of 
new dust species. The 21 /im emission would be sup- 
pressed in regions where the smallest grains have been 
destroyed by the central star, and in regions of the dust 
shell where the temperature and density are too low. In 
this case, the best candidates for the production of the 
emission feature would be cool (~ 100 K), micron sized 
/3-Si C and nano-SiC g rains that contain carbon impuri- 
ties (|Speck et al.ll2005T and references therein). Here, the 
a-SiC and /3-SiC refer to different polymorphs of silicon 
carbide grains. 

The scenario outlined in ISpeck et alj (|2005l ) appears 
to be consistent with our observations of the IR shell in 
G54.1+0.3. The 21 /im feature in G54.1+0.3 is strongly 
pronounced in the densest region, i.e. the IR knot, while 
it is suppressed in the diffuse shell emission that has a 
lower density. Our observations seem to suggest that 
SiC grains are strong candidates for producing the ob- 
se rved emission featur es. However, based on the models 
of iKozasa et al.l (|2009f) that predict relative abundances 
of grain species formed in SN ejecta, Si02 is also a likely 
candidate. We cannot currently rule out either of these 



grain species. In order to identify the composition that 
produces the 21 /im feature in G54. 1+0.3, more detailed 
modeling of the IR spectrum is required. The dust in 
G54. 1+0.3 is likely freshly-formed dust in the SN ejecta, 
with a similar composition as in Cas A. 

5.2. Extended IR Knot 

The bright IR knot in the western part of the IR shell 
is shown in Figure \5jp, and in Figure [3] where the 24 
/im emission is shown in blue. The emission across the 
knot is not uniform, but appears to be concentrated in 
three separate clumps of emission whose spatial profiles 
are broader than the MIPS 24 /im PSF. The extended 
clumps do not have counterparts at any of the IRAC 
wavelengths, nor does their emission correlate with the 
peak brightness of the MIPS 70 /im image. The high- 
resolution IRS spectrum of the knot is dominated by a 
broad emission feature around 21 /im and a rising under- 
lying continuum (Figure [8]) . The plot in Figure [TT] that 
depicts the spatial variation of emission features across 
the LL slit shows that the emission from the 21 /im fea- 
ture is an order of magnitude brighter at the IR knot 
than in the rest of the shell. We also note that the 9 and 
12 /im features are enhanced at the position of one of the 
point sources in the SL spectrum (Figure [13]) . If these 
features correlate with the 21 /im feature, then the 21 
/im feature may also be enhanced in the vicinity of the 
11 point sources due to a higher temperature, thus caus- 
ing the unusual [8] -[2 4] /i m color excess in the color-color 
plots of lKoo et all l|2008l ). 

Since photoionization by the PWN cannot adequately 
heat the dust in the IR knot, we explore the possibil- 
ity that this dust is shock-heated. The pulsar's X-ray 
jet appears to terminate at the location of the knot (see 
Figure [3]) , which may be an indication that it is driv- 
ing a shock into the IR shell and compressing the dust 
to produce a small dense region with enhanced 21 /im 
emission. The [S III] 18.7 /im emission and the ratio of 
the [S III] 18.7 /xm to 33.5 /im lines show a peak along 
the knot, implying a significant density enhancement in 
this region. The density calculated from the line ratios in 
the LH and SH at this position is in the 500-1300 cm~ 3 
range. In order to check if this scenario is plausible, we 
estimated the fraction of the pulsar's spin-down luminos- 
ity (E) that would have to power the jet in order to get 
enough ram pressure to compress the dust in the shell. 
The spin-down luminosity E for the pulsar in G54. 1+0.3 
is 1.2 x 10 37 erg s _1 . Assuming a a shock velocity of 25 
kms -1 fSection l4.5p . a preshock density of 30 cm -3 , and 
a mean particle weight of 30 m p , the shock dissipates 1/2 
pv 3 = 0.012 erg cm" 2 s — 1 . By estimating the area of 
the IR knot to be 600 square arcseconds (0.52 pc 2 at a 
distance of 6 kpc), we calculate that the power in the jet 
would need to be ~ 6 x 10 erg s _1 , or approximately 
half a percent of E, to produce the required ram pressure. 
This represents a reasonable (if approximate) fraction of 
the available spin-down power, indicating that the jet 
is compressing the dust in the shell and producing the 
denser IR knot with the pronounced 21 /im feature. 

5.3. Dust Mass Estimate 

Since the spectrum of the IR shell is dominated by the 
emission from the 21 fim feature, it was difficult to obtain 
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Fig. 12. — Intensity of the [Ar II], [Ne II], and 24 /im shell emis- 
sion as a fu nction of position across the SL slit. The slit is shown 
in Fig. 1111 Positions 0—30 correspond to positions across the SL 
slit from northern to southern edge of the IR shell emission. The 
statistical 1-tr uncertainties on each data point are on the order of 
15%. The pixel scale for the SL slit is 1.8" /pixel. 



a satisfactory fit to the spectrum with any single grain 
composition. In order to make a rough estimate of the 
dust mass in the IR shell, we fit the spectrum of the un- 
derlying continuum with two d ifferent grain compo sition 
models; as tronomical silicate s ([Laor fc Drainelfl993f ) and 
Mg 2 Si04 (|Jager et al.l 12003V We excluded the spectral 
region from 18-25 pm where the emission from the 21 
/im feature dominates. In the fitting, we used the high- 
resolution spectrum from the diffuse shell emission, posi- 
tion 1, where the 21 pm feature is somewhat suppressed, 
as well as the MIPS 70 pm data to help constrain the dust 
temperature. Before fitting, the spectrum was scaled so 
that the integrated flux over the MIPS 24 pm bandpass 
is equal to the total extinction-corrected flux at 24 /im of 
40 Jy. The fit was not satisfactory since it left large resid- 
uals from the 12 and 21 /im features, but it can still be 
used to make a rough estimate of the amount of silicates 
that are present in the IR shell. Grains of astronomical 
silicates at a temperature of ~ 68 K, and Mg2SiC>4 at 
~ 63 K, can approximately match the continuum and 
the 70 /im data point, excluding the 18-25 /im region 
and broad emission features at shorter wavelengths. The 
total dust mass in the shell is given by 



Mdust = 



F u d 2 Apa 

B v (Td) ZQabs 



(1) 



where F„ is the total IR flux, d is the distance, B„ 
is the Planck function evaluated at the grain temper- 
ature, p and a are grain density and size, and Q a bs 
is the absorption efficiency. Absorption efficiencies for 
0.05 am sized astro nomical silicates and forsterite are 
from lLaor fc Draind (|1993h and Uager et~aTI (|2003h . re- 
spectively. The assumed distance is 6 kpc. We find an 
approximate dust mass of 0.02 M Q for astronomical sil- 
icates, and 0.04 M Q for forsterite. Since the fits were 
not satisfactory for all of the observed emission features, 
these values should only serve as rough estimates. We 
note that this estimate is in the sa me range as the dus t 
mass estimated for the Crab Nebula (jTemim et al.l f2006). 




10.5 11.5 
Microns 

Fig. 13. — Sample spectra from three different positions along 
the SL slit. The spectra were extracted from 6 pixel long sub-slits 
at positions showing strong and weak [Ne II] emission, and at the 
position coincident with an embedded point source, where the 9 
and 12 (im broad emission features are most prominent. The pixel 
scale for the SL slit is 1.8"/pixcl. 

6. PHYSICAL INTERPRETATION FOR THE IR SHELL 

The cartoon in Figure [14] summarizes the IRS spec- 
troscopic results for different regions of the IR shell sur- 
rounding the PWN in G54.1+0.3. The PWN fills the 
cavity of the IR shell, with the diffuse X-ray emission 
extending almost to the outer boundary of the MIPS 24 
/im emission (Figure [3]) . The dark grey region in Figure 
1141 represents the location where we observe enhanced Si 
and Ar abundances, and where the intensity of the silicon 
line shows a sharp rise (see Figure [TT] ) The highest line 
broadening is also observed in this region (Table [5J . The 
diffuse shell emission, represented by light grey, shows 
lines that are somewhat less broadened and continuum 
emission with a weak 21 /im dust feature. The pulsar's 
jet appears to terminate at the bright IR knot whose 
spectrum shows strong emission from the 21 /im feature 
and a higher density, as derived from the S III 18.7/33.5 
pm ratio. 

In this section, we consider two possible physical sce- 
narios that may explain the observations of the IR 
shell. A previous explanation is that the shell formed 
by the progenitor's wind that compressed the surround- 
ing medium and then fragmented and c ollapsed to form 
YSOs, as suggested bv lKoo et al.l (|2008| ). An alternative 
scenario is that the shell is composed of expanding SN 
ejecta, and that the observed point sources are produced 
by radiative heating of ejecta dust by early-type stars 
embedded within the expanding SNR. Based on the IR 
data, and the lack of thermal X-ray emission from the 
shell, we find the ejecta interpretation a more plausible 
explanation for the origin of thermal dust emission in 
G54.1+0.3. 



6.1. Preexisting Shell 
The 24 pm MIPS image of the IR shell in G54.1- 



0.3 



contains 11 point sources embedded in the diffuse shell 
emi ssion that a r e arra nged in a ring-like structure (Figure 
S]). iKoo et ail (|2008| ) suggest that these point sources 
are massive pre-main-sequence stars, larger than 10 Mq, 
with an age of < 2 Myr. Their formation is suggested to 
be triggered by the progenitor star of G54. 1+0.3 during 
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its post-main-sequence stage of evolution. Most of the 
point sources are detected in IRAC and in the near-IR 
2MASS images, and some even have optical counterparts. 
Their JHK S color-color and color-magnitude diagrams 
indicate that they are young stars with spectral types 
r anging from Bl. 5 to 08. 

iKoo et a l. (2008) show that the spectral energy distri- 
bution (SED) of the brightest star in the K s band, the 
southernmost point source in Figure [H is similar to the 
spectral energy distribution of early Herbig Be stars. The 
IRAC and MIPS 24 /im colors, however, are not consis- 
tent with the usual colors of Class 0/1 and Class II YSOs. 
While the stars in the IR loop have a small IRAC 3.6] 
- [5.8] excess, consistent with Class II sources, the [8] - 
[24] is much higher than exp ected, in the 6 to 8.5 magni- 
tude range (|Koo et al.l l2008) . The IR spectrum in Figure 
[8] shows that there is a broad emission feature peaking 
around 21 /im. This feature likely dominates the MIPS 
24 /im flux throughout the diffuse shell, and it may be 
the cause of the unusually high [8.0] - [24] excess. It may 
be that the point source photometry is contaminated by 
the emission from the diffuse shell, but since the observed 
24 /im excess is up to 4 magnitudes higher than what is 
expected for Class II YSOs, it appears more likely that 
the emission from the 21 /im feature is enhanced in the 
vicinity of the these sources. A similar enhancement is 
observed in the IR knot, discussed in Section [5T2l 

If the point sources in the IR shell are embedded YSOs, 
their presence would require that the shell was formed 
prior to the SN explosion. However, the high expansion 
velocity of over 500 km s" 1 , as derived from the silicon, 
iron, and chlorine lines, indicates that these lines have 
to come from rapidly expanding SN ejecta. Another in- 
consistency with the interpretation of a preexisting shell 
is the absence of thermal X-ray emission in G54. 1+0.3, 
and a lack of evidence for any interaction with the super- 
nova blast wave. If the SNR has encountered the shell, 
we would expect the reverse shock to have disrupted the 
PWN. The structure of the PWN shows no indication 
that the reverse shock has reached its surface, providing 
further evidence against this scenario. In the following 
section, we provide an alternative interpretation for the 
origin of the IR point sources and the diffuse shell emis- 
sion, without the need for a preexisting shell. 

6.2. SN Ejecta 

Here, we consider a scenario in which the IR shell is 
composed entirely of SN ejecta. The shock models im- 
ply that the PWN overtakes the freely-expanding ejecta 
with a shock velocity of ~ 25 kms^ 1 (see Section |4~5")) . 
The line emission from the IR knot, where the pulsar's 
jet compresses the ejecta, is likely dominated by shock 
emission. However, in the rest of the shell where we only 
have upper limits on the density, photoionization by the 
stellar radiation field and the PWN likely accounts for 
most of the emission. If the density of the shell is signif- 
icantly lower than in the IR knot, the surface brightness 
of [Si II] produced by slow shocks would be too low to 
match the observed [Si II] intensity. However, if the den- 
sity of the shell is comparable to the density of the IR 
knot, the ram pressure of the PWN would be too low 
to drive a 25 km/s shock into such dense ejecta. Emis- 
sion from the shock most likely dominates only in the 
region where the jet drives a slow shock into the dense 
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Fig. 14. — Summary of the IRS spectroscopic results for the IR 
shell in G54.1+0.3. 

IR knot region, while the emission from photoionization 
dominates in the rest of the shell. The sharp [Si II] peak 
at the inner edge of the shell (Figure fTTj) may originate 
from a region with a density that is lower than in the 
IR knot, but high relative to the rest of the shell. This 
region may be produced by the compression of gas by 
the PWN, but with a lower ram pressure than in the jet. 

The high inferred preshock ejecta density in the IR 
knot implies that dense ejecta clumps with a low vol- 
ume filling fraction are embedded in less dense in- 
terclump medium. Such clumps are expected to be 
produced by compression of ejecta by a radioactively- 
heated and expanding Ni bubble. A s imilar scenario 
has b een suggested for SNR 0540-69.3 (|Williams et all 
120081 ). From equation (5) in Williams et al. (2008) 
that is valid for a Type IIP ejecta model, the pre- 
dicted densities in a shell swept up by the Ni bubble are 

- 250(t S NR/1500yr)- 3 (M j/15M Q ) 5 / 2 E 51 3/2 amucm" 3 . 
At 1500 yr and with M ej = 2OM , E 51 = 1, this gives 
500 amu cm -3 . Assuming a mean particle mass of ~ 30 
m p , and a preshock clump density of ~ 30 cm -3 , this is 
in reasonable agreement with observations. 

The broadened line emission suggests an ejecta veloc- 
ity greater than ~ 500 kms , since we are observing 
regions away from the center of the remnant where the 
expansion is mostly transverse to the line of sight. The 
point sources observed at 24 /im can be explained if we 
consider a scenario in which the SN exploded inside a 
cluster of young stars. Assuming a distance of 6 kpc, the 
observed stellar density in the region co vering the IR shell 
is no t unusual for a young cluster (e.g. lElmegreen et al.l 
2000). In this scenario, the SNR blast wave has expanded 
beyond some of the cluster stars, embedding them within 
the expanding SN ejecta. The ejecta dust is then blow- 
ing by the embedded stars and being radiatively heated. 
The dust model described in Section 16.2.21 shows that if 
the dust is heated by a dozen such stars, this can repro- 
duce the observed 24/70 /im ratio and give rise to the 
IR point sources at 24 /im without a need for triggered 
star formation. The fact that the spectrum of the shell 
closely resembles the sp ectrum of freshly-formed dust in 
Cas A (|Rho et al.l 120081 ) provides additional support for 
this interpretation. 
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6.2.1. Origin of Multiple Emission Regions 

As discussed in the previous section, the line emission 
from the shell likely arises from a combination of slow 
shocks driven into dense ejecta and photoionization. The 
IR spectra do indeed show that different emission lines 
have different spatial and kinematic properties (Section 
14. 2| and that they may arise from distinct components 
in the shell (Section 14. 3(1 . Intense UV stellar radiation 
within the stellar cluster may be heating the ejecta gas 
either by directly ionizing it or by ejecting photoelectrons 
from grains. Ionization of H and O requires photons with 
energies greater than 13.6 eV, that may be scarce if only 
B-type stars were present within the cluster. Abundant, 
less energetic UV photons that heat the dust also eject 
photoelectrons. Ejecta gas heated by grain photoelec- 
trons is expected to cool by emission in IR lines of singly 
ionized species with ionization potentials less than 13.6 
eV, including [Si II] 34.8 //m, [Fe II] 26.0 /im, and [CI II] 
14.4 /im. Unlike Si and Fe, CI is not expected to be 
strongly depleted onto dust, and this minor ejecta con- 
stituent may be the major coolant in freely-expanding 
ejecta. 

Heating by stellar photons should be effective both in 
the ejecta shell swept up by the PWN and in the freely 
expanding ejecta ahead of the shell. At sufficiently large 
radial distances, free-expansion velocity must exceed the 
shell velocity. If heating were still effective there, widths 
of [Si II] , [Fe II] , and [CI II] lines could be be larger than 
widths of lines such as [S III] that are more efficiently 
produced by shocks in the swept-up PWN shell. This 
is what is observed (Table [4|), suggesting that heating of 
ejecta gas by grain photoelectrons, and perhaps also by 
direct photoionization by more energetic photons, occurs 
in G54.1+0.3. 

6.2.2. Origin of IR Point Sources 

Here, we describe a model in which the ejecta dust is 
being heated by stellar sources inside a cluster in which 
the SN exploded. The stellar radiation field decreases 
in intensity as r~ 2 with increasing distance r from a 
star (and steeper still if the dust optical depth r is non- 
negligible), so dust is heated to much higher tempera- 
tures close to the star than far away from it. Hotter 
dust re-radiates the absorbed stellar radiation at shorter 
wavelengths than the more distant cooler dust, an effect 
that produces an apparent infrared excess in the stellar 
spectrum. We demonstrate this effect with a simple dust 
model, which is also used to estimate the total ejecta 
dust mass in G54. 1+0.3. 

We consider a BOV star with a luminosity of 25,000 Lq 
and temperature of 30,000 K, surrounded by uniformly 
distributed ejecta dust. There is a central cavity devoid 
of dust in the immediate vicinity of the star, where dust 
is cleared eithe r by a stellar wind or by the st ellar radi- 
ation pressure (|Artvmowicz fe Clampinl 119971 ). but this 
cavity is expected to be small in view of the large ram 
pressure exerted by fast-moving SN ejecta. We assume 
a small cavity radius of 0.003 pc; properties of the mid- 
IR emission are insensitive to the assumed cavity radius 
as long as the cavity remains small. Absorption of stel- 
lar UV radiation by dust and its re-radiation as ther- 
ma l dust emission in the IR are modeled as described 
by iBorkowski et all (|1994f) . These dust models include 
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Fig. 15. — Normalized model profiles in the 24 /im and 70 /im 
MIPS bands (left and right solid lines, respectively). The corre- 
sponding Spitzer's PSFs are plotted by dashed lines. Because of a 
large radial gradient in dust temperature, IR emission is spatially 
extended at 70 (im, but a point-like emission enhancement appears 
at 24 /im and at shorter IRAC wavelengths. 

effects of grain temperature fluctuations that are im- 
portant for small grains heated by energetic UV pho- 
tons. Scattering of photons by dust is not included, 
so these models underestimate absorption if t is not 
small. Mu ltiple dust species are predicted to form in 
SN ejecta ()Kozasa et all 120091) , one of the most abun- 
dant being forsterite (Mg 2 Si04). We considered amor- 
phous forsterite grains in our mod els, with density of 3. 3 
g cm" 3 and optical constants from lScott fe Dulevl (|1996j) . 
As a single grain size is unlikely, a power law distribu- 
tion in grain sizes was assumed, with a power-law index 
of —3.5, ranging in radii from 0.001 /im to 0.25 /im. We 
chose 0.007 M Q pc" 3 as the spaced-averaged grain mass 
density; within a factor of a few, this choice is compati- 
ble with the dust mass estimate presented in Section 15.31 
and the volume occupied by the IR shell. 

Our dust model matches the extinction-corrected 
70/24 /im MIPS flux ratio of 1.9 if we choose the ra- 
dius to which the dust extends to be 0.69 pc. The 70/24 
MIPS flux ratio decreases outward, tracing a steep ra- 
dial dust temperature gradient, so models with a smaller 
(larger) radius under(over)predict this ratio. We con- 
sider only dust within a radius of 0.69 pc in our discus- 
sion of spatial IR profiles and of dust mass estimates. 
Within this radius, the infrared dust luminosity is 1200 
Lq, comprising only 4.8% of the stellar luminosity, so 
the SN ejecta are optically thin to stellar radiation. The 
predicted surface brightness in the 24 /im and 70 /im 
MIPS bands is shown in Figure [T51 without accounting 
for the limited spatial resolution of Spitzer. A sharp, 
point-like emission enhancement is present at 24 /im, 
and absent at 70 /im. A comparison with the Spitzer 
PSF (Figure [T5|) reveals that the 24 /im emission en- 
hancement at the center cannot be spatially resolved with 
Spitzer. The uniformly-distributed ejecta dust appears 
as an unresolved point source, surrounded by faint ex- 
tended emission. This is no longer true at 70 /im; emis- 
sion at this wavelength is predominantly produced far 
from the star, and Spitzer's spatial resolution becomes 
adequate for mapping the spatially-extended dust emis- 
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sion. This simple dust model matches well the observed 
morphology of G54. 1+0.3 as seen by the IRAC, MIPS 
and Akari detectors (|Koo et alll2008l ). 

The 70 /im flux, equal to 7.5 Jy in our model with a 
radius of 0.69 pc, is produced by 0.0097 Mq of forsterite 
dust. Since the total spatially- integrated 70 /im flux is 
equal to 76 Jy, the 11 stars w it h an apparent infrared 
excess identified by iKoo et all (|2008l ) are sufficient to 
heat ejecta dust with an estimated total mass of approx- 
imately 0.1 M & . This mass is within about a factor of 
two of the mass estimate from forsterite grains, derived 
from spectral fitting is Section 15.31 This dust mass esti- 
mate can be improved once we learn more about these 
stars from their near IR spectroscopy, and once the dust 
composition is better constrained. In particular, a better 
understanding of the stellar radiation field in the parent 
stellar cluster (or association) of the G54. 1+0.3 progeni- 
tor will allow for more realistic modeling of dust heating 
by stellar UV photons than currently possible. 

7. CONCLUSIONS 

In this paper, we presented deep Chandra observations 
of G54. 1+0.3 and Spitzer imaging and spectroscopy of 
the surrounding IR shell. The 300 ks Chandra observa- 
tion was used to derive a new value for the absorbing 
column density toward G54. 1+0.3 that was then used 
in the IR analysis. A more detailed X-ray study will 
be presented in a separate publication. We fitted the 
X-ray spectra from different regions of G54. 1+0.3 with 
an absorbed power-law model and found that a more 
accurate value of N# results when the piled-up pulsar 
region is excluded from the fit. The new value of Njj is 
(1.95+0.04) x 10 22 cm -2 , somewhat larger than the previ- 
ous estimate by lLuet all (|2002t ) of (1.6±0.1)xl0 22 cm" 2 . 

The IRAC and MIPS images reveal an IR shell with a 
dozen point sources arranged in a ring-like structure, and 
a dense extended region that is aligned with the pulsar's 
jet. IRS spectra reveal a number of ionic emission lines, 
including [Ar II], [Ne II], [CI II], [S III], [Fe II], and [Si II], 
and a rising continuum with a broad emission feature 
at 21 /im. We conclude that the emission from the IR 
shell likely arises from SN ejecta and freshly-formed SN 
dust. The evidence for this interpretation comes from 
the observed morphology of the PWN and the shell, sig- 
nificant broadening of IR emission lines, comparisons of 
line intensities with shock models, and the observed fea- 
tures in the dust emission spectrum. Since the PWN 
fits into the cavity of the IR shell, these components are 
almost certainly associated and likely interacting. The 
spectral lines show evidence for broadening with a maxi- 
mum FWHM of over 1000 kms -1 , suggesting that these 
lines arise in rapidly expanding SN ejecta. Based on our 
shock model results and density diagnostics, the observed 
line intensities are consistent with a scenario in which the 
emission from the IR knot is produced by the pulsar's jet 
driving a ~ 25 kms" 1 shock into a dense ejecta clump. 
The emission from the rest of the shell is most likely dom- 
inated by photoionization of the ejecta gas. The [S III] 
18.7/33.5 /zm line is consistent with a low density limit 
in the IR shell and a higher density in the IR knot of ~ 



1000 cm -3 . This extended IR knot may be produced as 
the pulsar's jet compresses the shell material. In order 
for the jet to have enough ram pressure to produce the 
density enhancement, it would need to be powered by 
approximately 0.5 % of the pulsar's spin-down luminos- 
ity E. If the density in the rest of the IR shell is near the 
estimated upper limit, i.e. comparable to the density of 
the IR knot, then the shell would have to be composed 
of dense ejecta clumps with a low volume filling fraction, 
as produced by a radioactively heated and expanding Ni 
bubble. 

The broad 21 /im emission feature is enhanced at the 
position of the knot and weak in the rest of the shell. 
The profile of the f eature is remark ably similar to the 
spectrum of Cas A (Rho et al. 2008), and is most likely 
produced by the same dust species. T he emission may b e 
caused by Si02 grains, as sugg ested bv lRho et all (f2009h . 
but SiC grains are also likely candidates since they are 
capable of producing the 21 /im feature and other associ- 
ated emission features around 9 and 12 /im ([Speck et all 
2005). We estimate the dust mass in the IR shell to be in 
the 0.04-0.1 Mq range, assuming a forsterite grain com- 
position. More detailed modeling of the IR spectrum is 
required to better characterize the properties and com- 
position of the dust. 

We find that our observations are best described by a 
scenario in which the shell emission arises entirely from 
SN ejecta. In this case, the point sources in the shell 
are not young stellar objects, but would be attributed 
to radiatively heated ejecta dust by early-type stars that 
belong to a stellar cluster in which the SN exploded. A 
simple dust model shows that this scenario can reproduce 
the observed IR emission, including the IR point sources. 
The stars are illuminating unshocked dust grains and we 
are probing freshly formed SN dust before its encounter 
with the reverse shock. The study of dust production in 
SNe may prove to be most effective in cases where the 
SN explosion occurs inside a young stellar cluster. 
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